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Early structure formation in quintessence models and its 
implications for cosmic reionisation from first stars 
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ABSTRACT 

We present the first hydrodynamic N-body simulations of primordial gas clouds responsible 
for the reionisation process in dark energy cosmologies. We compare the cosmological con- 
stant scenario with a SUGRA quintessence model with marked dynamics in order to highlight 
effects due to the different acceleration histories imposed by the dark energy. We show that 
both the number density of gas clouds and their dumpiness keep a record of the expansion rate 
during evolution, similar to the non-linear dark matter profile at virialisation, as was recently 
demonstrated by Dolag et al. (2004). Varying the shape of the primordial power spectrum, we 
show how this effect is mitigated by a running spectral index decreasing the power at small 
scales. Our results demonstrate that, in order to constrain the dark energy from large scale 
structures, one must track its effects down to the distribution of luminous matter. 

Key words: early universe - cosmology: theory - galaxies: formation 



1 INTRODUCTION 

Different observational dat a, like those fr om high-redshift super- 
novae jRiess et all2004l Astier et al. 2006), the cosmic microwave 
background l Spergel et al. 2003, 2006), and large-scale structure 
I Tegmark et al. 2004 ; Cole et al. 2005), are now giving a consistent 
picture of our universe, which can be described by the so-called 
'concordance' ACDM model: a spatially flat universe whose ex- 
pansion accelerates in the present epoch because of a dominant 
dark energy component. Historically, the first and simplest candi- 
date for that has been the cosmological constant A, correspond- 
ing to matter with an equation of state p — wpc 2 , with constant 
W = — 1. However, observations suggest a value for A which is 
more than a hundred orders of magnitude smaller than the energy 
scales expected to be responsible in the very early universe. 

To avoid inelegant solutions based on parameter fine-tuning, 
the general idea of dark energy is extended to encompass the 
so-called quintessence, possibly corresponding to a suitably self- 
interacting scalar field, whose pressure and energy density evolve 
during cosmic history. The currently available observational data 
sets do not allow yet to place strong constraints on w. We know that 
it has to come close to — 1 with a precision o f about 10 per cent a t 
the present epoch (see, e.g. jRiess et all2004HSpergel et alj 2006). 



but its redshift evolution is still poorly constrained and can hope- 
fully b e determined only with new-generation data; see Seli ak et alJ 
1 2005) for one of the very first attempts to measure the high-redshift 
dark energy equation of state. 

One consequence of tu > —1 earlier is that the forma- 
tion of cosmic structures sets in earlier compared to the ACDM 
model. At a given redshift, this corresponds to a hi gher abun- 
dance of more concentrated halos (see, e.g.. iDolag et aljl2004l) . 
This may open alternative ways to stu dy quintessence models based 
on galaxy-cluster counts (Vfa n fi ^ Steinhardt 1998; Haim an et alJ 
1200 ltlBattve & Welle j2003UMaiumdar & MohJ2003l) and strong 
gravit ational lensing l Bartelma nn et alJ 120031 : Meneghet ti et alJ 
2005). 

The different history of structure formation must also affect 
the reionisation epoch. To produce the required ionising radiation, 
a sufficiently high number of early stellar sources is needed, such 
as super-massive Pop-Ill stars or more "standard", massive Pop-II 
stars in proto-galaxies. In particular, Pop-Ill stars are supposed to 
be composed mainly of hydrogen and helium with their primordial 
abundances, to have masses > 100 Mq much larger than the stan- 
dard Pop-I and -II stars, and to form in dark matter haloes with 
typical masses of 10 5 - - 6 M Q (the so-called mini-haloes). 



2 U. Maio et al. 



The recent analysis of the temperature-polarisation cross- 
correlation and the polarisation auto-correlation measured by the 
WMAP satellite now allows the estimation of the epoch when 
reionisation occurs. The Thomson optical depth of of r ~ 0.17 ± 
0.04 extracted from the first-year WMAP data had been inter- 
preted as signalling the high redshift of z = 17 ± 5 fo r the pri- 
mordial epoch of global reionisation iKoeut et aljl2003h . The re- 
cently released three-year WMAP data allowed an improved con- 
trol in particular of the polarised foreground emission, yielding 
r = 0.09 ± 0.03. This corresponds to z = 10^ '3 f° r me comple- 
tion of the reionisation process, even if some level of parameter de- 



ferent authors (see, e.g., 


Wvithe & Loeb 2003; Ciardi et al. 2003; 


Sokasian et al. 2003, 


2004 


1 suggested that these reionization data 



can place new complementary constraints on the cosmological pa- 
rameters, and in particular on the nature of dark energy. 

So far, studies of the structure formation process in dark 
energy cosmologies were focused on modifications of dark mat- 
ter structures due to the different expansion histories (see, e.g., 
iKlvpin et alj2003l:lDolag et al.l2 004V They revealed that virialised 
objects keep a record of the expansion rate at the time of their for- 
mation. The general picture is that in cosmologies with w > — 1, 
the increase of the dark energy density with redshift enhances struc- 
ture growth, and virialised haloes become mo re conce ntrated be- 
cause of the denser environment they form in iDolag et alJl2004|) . 
This effect is stronger in tracking quintessence models, which we 
describe in the next section, in which the increase of the dark energy 
density with redshift is enhanced. Complementary hydrodynamical 
simulations still need to be carried out. This work is a first step into 
this direction, beginning with the numerical study of the depen- 
dence of the primeval gas clouds responsible for the reionisation 
process on the cosmic expansion rate at their formation time. 

Due to the higher concentration of structures in dark en- 
ergy scenarios, which has so far been verified only for the dark 
matter, it is mandatory for our study to control the shape of the 
primordial power spectrum. When combined with the Ly-a for- 
est data and with the analysis of the 2dF galaxy redshift survey, 
the first-year WMAP results supported a "running" spectral index 
which tilts the spectrum slightly towards small scales, starting at 
k > 1 Mpc" 1 . The thre e-years WMAP data are compatible with 
this JSpergel et alj200r3) . albeit with less emphasis than the earlier 
results. A running index would cause a slower growth and evolu- 
tion of small-scale cosmic structures compared to models with a 
constant index n. A reduction of power on small scales may al- 
leviate several potential discrepancies in the ACDM model, such 
as the abundance of substructures in galaxies and the high central 
concentration of galactic hal oes. However, as shown by an alytic 
and n umerical work (see, e.g. lSomerville et alj2003llYoshida et alJ 
003b), models with a running spectral index (RSI) may have severe 
problems in producing enough objects to allow a global reionisa- 
tion at high redshift. Since this might be balanced by the enhanced 
structure growth due to the dark energy, it is important to jointly 
study these two aspects. 

In the present paper, we study the high-redshift structure 
formation in quintessence models based on the results of high- 
resolution cosmological N-body simulations combined with hydro- 



1 We have also to notice that the three-year WMAP data are suggesting a 
lower value not only for t, but also for the power spectrum normalization 
erg: these two effects nearly cancel in terms of early structure formation, i.e. 
the inferred small r does not really allow slow reionization for a given erg . 



dynamics. Specifically, we consider four different cases combining 
two flat cosmological models, i.e. the concordance ACDM model 
and a quintessence model with a SUGRA potential, with two types 
of the primordial power spectrum, one with a constant spectral in- 
de x n = 1 and a n RSI model assuming the best fit relation found 
bv lSpergel et alJ 1 2003). More detail will be given in the next sec- 
tion. The paper is organised as follows. In Sect. [2] we describe the 
general characteristics of the hydrodynamical simulations used be- 
low and introduce the quintessence models chosen. The techniques 
adopted for identifying dark matter haloes and the corresponding 
results are presented in Sect. [3] The abundances of gas clouds are 
presented in Sect.|4]together with the clumping factors and recom- 
bination times as computed from the simulation outputs. The impli- 
cations of the previous results in terms of reionisation are discussed 
in Sect.|5| The final discussion and our main conclusions are drawn 
in Sect.© 



2 THE SIMULATIONS 

We first review the cosmologies chosen, and then the computa- 
tional and hydrodynamical features of our numerical simulations. 
As mentioned in the introduction, this work is a first step towards 
investigating the dependence of hydrodynamic N-body simulated 
structures on the global cosmic expansion rate, and it is appropri- 
ate to study the connection with the shape of the primordial power 
spectrum. Thus, we need dark energy cosmologies with different 
dynamics as well as constant and running index in the primordial 
power spectrum. 

The standard, concordance ACDM model is our reference 
case. It is spatially flat, has the present matter density f2om = 0.3, 
the baryon density J7ob = 0.04, and a dominating cosmologi- 
cal constant, 57qa = 0.7. The Hubble constant is assumed to be 
h — 0.7 in units of 100 kms" 1 Mpc~\ The second model is cho- 
sen so as to emphasise the dynamical effects of the dark energy. 
The spatial curvature is also vanishing, with f2om = 0.3 and 70 
per cent of the critical density contributed by a quintessence scalar 
field $ evolving in time under the effect of a potential energy V(&). 
We review here t he main characteristics of this model, referring to 
IPeebles & Ratrl 1 2003 ) and references therein for further detail. 

The field $ must satisfy the Klein-Gordon equation 



$ + 3H$ + ^ 



: 



The equation-of-state parameter w = pj (pc 2 ) is 



2^ 



i$ 2 + T/($) 



(1) 



(2) 



Evidently, w evolves with redshift, but the behaviour of the cos- 
mological constant is reproduced as m — > —1 in the static limit 
$ 2 < V($). 

The attractive feature of such models are their attractor solu- 
tions which exist for exponential and inverse power law potentials 
or combinations of those. Suitably tuning the potential amplitude, 
they can reach the present level of dark energy starting from a large 
range of initial conditions. The evolution of the equation of state 
may be more or less pronounced along these trajectories. We wish 
to have pronounced dark-energy dynamics in order to clearly iden- 
tify its effects. Thus, we choose a SUGRA potential 

= — exp(4fG$ 2 ) , 



(3) 
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Figure 1. Upper panel: The redshift evolution of the equation-of-state pa- 
rameter w for the SUGRA model with wo = —0.85 (solid line), compared 
to the concordance ACDM model (dash-dotted line). Lower panel: The lin- 
ear growth factor (normalised to the ACDM model at present) for the same 
cosmologies as in the upper panel. 
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Figure 2. Comparison between the two different CDM power spectra 
adopted in generating the initial conditions for the numerical simulations: 
solid and dotted lines show spectra with fixed n = 1 and with the run- 
ning spectral index J4}, respectively. The arrow indicates the range of length 
scales covered by our simulations. 



where G is the Newton constant, as su ggested by super- gravity cor- 
rections iBrax & Martin 1999, 2000). By exploiting the properties 
of tracking trajectories, we fix M by requiring that the dark energy 
makes up 70 per cent of the critical density today, and a — 6.5 
so that w — > wo = w(z = 0) = —0.85 today. At higher red- 
shift, w > wo I Brax & Martin 2000), which meets our purpose of 
emphasising differences with the ACDM model. It should be noted 
that this model i s already at variance with t he observ ational con- 
strain ts (see,e.g., ISpergel et all2003l : lKnop et all2003tfRiess et alJ 
2004). 

The upper panel of Fig. \\\ shows the redshift evolution of w 
for this model. The differences with the standard ACDM model 
(constant w — —1) are large at all redshifts. In particular, dur- 
ing the epoch of reionisation which we are mostly interested in, 
the SUGRA model has w w —0.35. Notice that w strongly af- 
fects the formation history of cosmic structures through its mod- 
ified expansion history. This is evident from the lower panel in 
Fig.Q which shows the linear growth factor (normalised to coin- 
cide with the ACDM model at z — 0) for the model shown above. 
In quintessence models where the dark energy is more abundant in 
the past with respect to now, structures tend to form earlier if the 
perturbation power spectrum is normalised at the present epoch. 
This is due to the earlier dark energy dominance, slowing down 
the perturbation growth. In our SUGRA model with Wo = —0.85, 
fluctuations at z = 20 are expected to be ~ 15 per cent higher than 
in the ACDM model. Consequently, we expect that non-linearity to 
be reached earlier, with obvious consequences for the reionisation 
epoch. 

Dynamical dark energy generally incre ases the concentra- 
tion of virialised structures I Dolag et al. 2004), which can also be 
caused by a modified primordial power spectrum. Therefore, we 
vary the simulation set-up using two different shapes of the pri- 
mordial power spectrum. Specifically, for both cosmological mod- 
els (ACDM and SUGRA), we construct initial conditions corre- 



sponding to two scenarios, namely a power spectrum with the stan- 
dard constant n = 1 normalised to as = 0.9. The second model 
also has as = 0.9, but allows for a reduced primordial power on 
small scales as suggested by the combined analysis of the first-year 
WMAP data, galaxy surveys and Ly-a data (see, e.g. lSpergel et alJ 
2003 ), within the range compatible with the three-year WMAP data 
ISpergel et alfeOOd) . 

The case for a running spectral index is interestingly contro- 
versial. Re-analyses of the Ly-q data (see, e.g. IViel et alJ[2004t 
ISeliak et all2005tlViel & Haehnell2006tlZaroubi et all200fj) lend 
much less support for a running spectral index than reported based 
on the first-year WMAP data and favour a constant primordial spec- 
tral in dex very close to unity (see, however, iDesiacaues & Nusserl 
2005). On the other hand, a running s pectral index is still viable 
and fits the three-year WMAP data well l Spergel et al. 2006), while 
the suppor t of a pure power law from large scale structures persist 
(see, e.g.. IViel etai]l2006l Iseliak et alj|2006t) . As mentioned be- 
fore, we choose our cosmological models so as to cover a broad 
range of phenomenologies rather than to fit the existing data well. 
Thus, we choose to adopt a primordial power spectrum given by 
P(k) oc k n{k \ with 



n(k) = n(fco) + 



1 dn / k 

2 dink 11 Vfc^ 



where k = 0.05 Mpc" 1 , n(Jfco) = 0.93 and 
dn 

-0.03 . 



dink 



(4) 



(5) 



This choice co nveniently allows a di rect comparison with the sim- 
ilar analysis bv lYoshida et alJ lo03bl) . who adopt the same primor- 
dial power spectrum, but only consider the concordance ACDM 
model. It is worth recalling that this represents the best fit to the 
first-year WMAP data and r emains compatible w ith the three-year 
data, dn/dlnfc < -0.06 ISpergel et af1l2003L EoOol) . Figure |2| 
shows the two primordial power spectra considered here. Evidently, 
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they differ at large k where the RSI spectrum is steeper. Less power 
on small scales will delay the structure formation (see also the cor- 
responding discussion in Yoshida et al. 003b). 

Combining two cosmological models (ACDM and SUGRA) 
with two primordial power spectra (constant n — 1 and RSI), we 
ran four different simulations; ACDM and SUGRA with n = 1 
and ACDM-RSI and SUGRA-RSI with a running spectral index. 
The initial conditions were generated at z ~ 120 using a comov- 
ing pe riodic box of 1 Mpc side length. We followed Dolag et al. 
1 2004) in adapting the initial conditions for ACDM to those for the 
SUGRA model. Thus, by construction, only the pairs of models 
sharing the same primordial power spectrum have the same ran- 
dom phases, allowing a direct comparison of the structures formed. 
The models with different primordial power spectra have different 
random phases and can thus not be directly compared (see Fig[3}. 

Since we are interested in structures during cosmic reionisa- 
tion, the box size used here can be considered sufficiently large to 
suppress cosmic varian ce (see the discussion on effe cts of finite box 
sizes in Yoshi da et al1l003bt iBasla & Prasadl f2006). Moreover, to 
prevent the non-linear regime from reaching the fundamental fluc- 
tuation mode in the simulations, we stopped the simulation at the 
latest at z ~ 15. 

When the gas starts cooling significantly within the first form- 
ing structure, the simulation time step can become very small (since 
we neglect feedback in these simulations), effectively stopping the 
simulation. Accordingly, we had to stop some simulations at a 
higher redshift, for example the SUGRA simulation, which could 
only be evolved until z ~ 19 because of its more rapid evolution, 
as we shall see in the later analysis. Nonetheless, all simulations 
are followed long enough for them to form the sufficient amount of 
structure we need for our analysis. 

The si mulations were performed using the GADGET2 
code fepringel et al] 1200 it ISpringellEool on the IBM-SP4/5 at 
CINECA (Bologna). GADGET2 is based on the combination of 
a tree-particle-mesh algorithm to solve the gravitational forces, 
and the smoothed particle hydrodynamics (SPH) scheme to de- 
scribe the hydrodynamics of gas particles. The code follows the 
non-equilibrium reactions of different chemical species (e~ , H, 
H + , H~, He, He + He ++ , H 2 , H? ), using the reaction coef- 
ficients computed by |Abel_etalJ ^997|) (for more detail, see also 
lAnninos et all Il997t; lYoshida etalT l003a). and adopts the cool- 
ing ra te due to molecular hydrogen estimated by iGalli & Pallal 
1 1998). The density field was sampled with 324 3 dark matter 
(DM) particles and an equal number of gas particles, having a 
mass of ~ 1040 Mq and 160 Mq, respectively. The comoving 
Plummer-equivalent gravitational softening length was fixed to 
e = 0.143 kpc. 

In Fig.|3|we show the projected gas and (mass-weighted) tem- 
perature distributions for all four models at z = 20.6. As ex- 
pected, the RSI simulations appear delayed, exhibiting smoother 
density and temperature distributions. The SUGRA models reveal 
their expected behaviour with a more evolved density field than 
their ACDM counterparts. Since the models with and without RSI 
have different random phases in their initial conditions, they are not 
supposed to reproduce the same structures. 



3 THE ABUNDANCE OF DARK MATTER HALOES 

The abundance of virialised dark matter haloes is usually expressed 
by the mass function N(M). From a theoretical point of view, 
different relations are available to predict the redshift evolution 



of N(M) once the cosmological model is fixed. Some of these 
were obtained from a specific model for the collapse of structures 
IPress & SchechteJll974l ISheth & Tormenll2002h : in other cases, 
they rep resent best fits to N-body simulations with very high res- 
olutio n iSheth & TormerJll999t |jenkins et ai]l200lk Iwarren et alJ 
2006). However, the theoretical mass functions have not yet been 
extensively tested in the regime of very low masses and very high 
redshifts which we are mainly interested in, where their predictions 
differ by orders of magnitude. This is due to the difficulty of com- 
bining sufficient resolution in space and mass with regions large 
enough to be considered a fair samples of our universe. 

The existing studies yield contradictory results. Simulations 
covering large volumes, such as those presented in Reed et al. 
(2003) (432 3 particles in a 50 h' 1 Mpc box) and in lSpringel et al 
1 20051) (2160 3 particle s in a 500 ft" 1 Mpc box), found that the 
I Sheth & Tormenll999L hereafter ST99) formula fairly reproduces 
the numerical results, but they can safely identify haloe s with 
masses > 10 10 M at z < 10. At z = 15, iReed et all <2003l) 
found a slight tendency of the ST99 model to over-predict the halo 
counts. 

Numerical work based on much smaller boxes, such as 
Ijang-Condell & Hemquisil 1200 J) (128 3 particles in a 1 h' 1 Mpc 
box) or l Yoshida et alJ J003bld> (2 x 324 3 particles in a 1 ft -1 Mpc 
box), found very goo d agreement with the predictions of the 
I Press & Schech teJl974 hereafter PS74) model. Notice that these 
simulations are more affected by their finite box size (see the 
following discussion). Further support for the PS74 formalism at 
high redshifts is provided by the analysis of iGao et alj _j200j). 
who s howed that the extended Press-Schechter theory iBoncLeHIl] 
1991) correctly predicts the structure growth and the dependence 
of halo abundances on the density of the surrounding region in the 
redshift range between 50 and 30. One of the best compromises 
between mass r esolution and volume coverage was achieved by 
iHeitmann et alJ 120061) . who used an extended set of simulations 
with 256 3 particles in boxes sized between (4 . . . 126) hT 1 Mpc. 
Their mass function at high redshift turns out to deviate signifi- 
cantly from the PS7 4 prediction for M > 10 7 h' 1 Mq, while the 
agreement with the iWarren et all200d. hereafter W06) relation is 
better. Notice that the W06 formula predicts a substantial suppres- 
sion of the halo abundance in the high-mass regime compared to 
the ST99 model. Finally, using an array of high-resolution N-body 
simulations with box sizes in the range 1 - 3000/i _1 Mpc, Reed 
et al. (2006) recently determined the mass function of dark matter 
haloes at redshifts between 10 and 30. They found that the PS74 
models gives a poor fit at all epochs, while the ST99 model gives a 
better match to the simulations, but still overpredicts the abundance 
of most massive objects by up to 50 per cent. 

We identify dark matter haloes in our simulations by running 
a friends-of-friends algorithm on the dark matter particles only and 
setting the linking length to 20 per cent of the mean inter-particle 
separation. The resulting mass of each halo is then corrected by the 
factor (1040 + 160)/1040 w 1.15 to account for the additional 
contribution by the gas particles. Figure|4]displays the redshift evo- 
lution of the number of haloes with masses exceeding 7 x 10 Mq 
for all four simulations. The figure compares the results with the 
predictions of the PS74, the ST99 and the W06 models, indicated 
by solid, dashed and dot-dashed lines, respectively. Note that we do 
not consider the relation found by Jenkins et al. (2001) because it 
cannot extrapolate to the mass and redshift ranges considered here. 

At first sight, the dark energy dynamics as well as the modi- 
fied initial power spectrum generate significant differences. First, as 
expected, the number of massive haloes is always lower in the RSI 
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Figure 3. The logarithm of the projected gas-density and mass-weighted temperature at z = 20.6 is shown in thin slices (l/7th of the box height) for the 
different models: SUGRA, ACDM, SUGRA-RSI and ACDM-RSI, as labelled in the maps. The density is in physical units and the slices are 6.6 (physical) kpc 
thick at z = 20.6. Note that the models with and without RSI have different random phases in their initial conditions and are thus not supposed to reproduce 
the same structures. 



models than in their scale-invariant counterparts. While the first ob- 
jects with M > 7 x 10 Mq start appearing near z ~ 30 in the 
SUGRA and the ACDM simulations, this is delayed to z ~ 25 and 
2 ~ 20 for SUGRA-RSI and ACDM-RSI, respectively. 

Second, the different expansion histories leave clear records 
in the structure population. Considering models with the same pri- 
mordial power spectrum, the halo abundance is always larger in the 
SUGRA than in the ACDM models by a factor ~ 2, confirming 
the earlier structure formation in dynamical quintessence models 
expected from the modified linear growth factor in Fig.Q 

We now consider the difference to the theoretical expecta- 
tions from the PS74, ST99 and W06 mass functions. All analytic 
relations were modified so as to account for the different expan- 
sion histories and initial power spectra in the different cosmologi- 
cal models. For the two simulations with a cosmological constant 
(ACDM and ACDM-RSI), the results indicate that the W06 for- 
mula agrees better with the simulations than the ST99 and PS74 
formulae. Some small deviations are seen at very high redshifts 
where, however, counts are low and the statistical uncertainty is 
high. The ST99 mass function tends to slightly overestimate the 
simulation results, again mainly at high redshift, where the dis- 
crepancies with W06 are more evident. On the contrary, the PS74 
mass function always underestimates the halo abundances, with 
differences reaching up to a factor ~ 3. Thus, our results for 
the cosmological- constant models agree with the previous analy- 
ses bylReed et alJj2003l) . ISpringel et alJ<2005ft and Heitman netai] 
(2006). However, this conclusion is in conflict with ot her studies 
Ijang-Condell & Hernauisl200lUYoshida et all003blci) . 

The situation is different for both SUGRA models. The mass 
function in the simulation with constant n = 1 is always larger 
than all theoretical predictions. In particular, the ST99 and W06 re- 
lations give very similar estimates which fall below the simulation 



results by 30 per cent, while the difference with PS74 is always a 
factor of 3-4. The tendency of PS74 to underestimate the numerical 
mass function is also confirmed for SUGRA-RSI (a factor of 3-5 
everywhere), while the agreement for W06 and mostly for ST99 is 
substantially better. The deviations reach 20 per cent only at low 
redshift. These results show that caution must be applied when ap- 
plying the standard theoretical mass functions to very high redshift 
in quintessence models. 

We must finally remark that, owing to the much higher mass 
resolution, the mass range covered in our simulations is quite differ- 
ent compared to earlier work. At the same time, our results could be 
affecte d by the small box size. In fact, as discussed in Yosh ida et alJ 
1 003b), the "effective" variance, i.e. the amount of power on the 
scales covered by the simulation, can be significantly smaller than 
that obtained by integrating the power over all scales. This dif- 
ference can systematically delay structure formation and possibly 
cause the discrepancies between the simulation results from differ- 
ent box sizes. To quantify the effect of the finite box, we apply the 
analytic expressions proposed by Baela & Prasad 1 2006), which al- 
low an evaluation of the required correction to the mass function. 
For example, assuming the PS74 approach and the redshift range 
covered by the data shown in Fig. [4] we find that the correction 
for M = 7 x 10 s Mq leads to a slight underestimate of 10-15 
per cent, which does not significantly affect our previous results. 
Similar conclusions can be reached on the ST99 relation. 



4 THE PROPERTIES OF THE GAS DISTRIBUTION 

We now show our results on the statistics and the physical proper- 
ties of the gas distribution. We concentrate on two aspects, namely 
the formation time of the gas clouds and their clumping factor, fo- 
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Figure 4. Number of haloes with mass exceeding 7 X 10 5 Mq per unit comoving volume as function of redshift in the different simulations (filled circles and 
dotted lines): SUGRA (upper left panel), ACDM (upper right panel), SUGRA-RSI (lower left panel), ACDM-RSI (lower right panel). Theoretical predictions 
based on the Press & Schechter (PS74, solid line), Sheth & Tormen (ST99, dashed line) and the Warren et al. (W06, dot-dashed line) relations are shown as 
well. 



cusing on their dependence on the underlying dark energy cosmol- 
ogy and primordial power. 

4.1 The formation of gas clouds 

We first study the formation of the gas clouds which are responsi- 
ble for reionisation. For this purpose, we run a friends-of-friends 
algorithm to group cold and dense gas particles, assuming a link- 
ing length of l/20th of the mean inter-particle separation to locat e 
the densest gas concentrations. Following Yosh ida et alj l003blct) . 
we only consider clouds exceeding 3000 M© , which corresponds 
to groups with at least 19 gas particles in our simulations. 



The results are shown in Fig. [5] in terms of the redshift evo- 
lution of the number of gas clouds n c iouds identified in the four 
simulations. It is not surprising that the trend between the models 
is similar as for the dark matter haloes in Fig.|4| Gas clouds start 
forming earlier and are more abundant in the SUGRA than in the 
ACDM model with the same primordial power spectrum. 

The suppression of power on small scales is, however, the 
dominant effect. In fact, the two RSI models have their structure 
formation delayed and have the smallest number of gas clouds in 
the final outputs. At z ~ 19, there are order-of-magnitude differ- 
ences between the models. We find ~ 300 clouds for SUGRA, 
~ 50 for ACDM, 7 for SUGRA-RSI and no clouds for ACDM- 
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Figure 5. Number of clouds per unit comoving volume as function of red- 
shift for SUGRA, ACDM, SUGRA-RSI and ACDM-RSI, as labelled. 



RSI. Also, the formation of the first gas clouds strongly depends 
on the cosmological parameters. It happens much before z — 30 
for SUGRA, while it is z ~ 18 for ACDM-RSI. The reason for 
this pronounced difference is the integral nature of the observable 
ficiouds considered here. The production rate of clouds per unit time 
is the key quantity, which differs between the cosmologies con- 
sidered because of the different expansion histories and primordial 
power spectra. The SUGRA model produces clouds earliest, fol- 
lowed by ACDM, while the same models with RSI stay in the same 
order but the structure formation process is delayed because of the 
modified primordial power. The integrated cloud count n c i ou( ja ac- 
cumulates these differences during the whole epoch of cloud pro- 
duction. 

Finally, we have to notice that in our study we neglected the 
feedback effects which might influence the structure formation his- 
tory of primordial clumps. For example, radiative feedbacks, such 
as photoionisation or molecule photo-dissociation, could stop cool- 
ing suppressing the formation of small objects; mechanical feed- 
backs and shocks can determine gas removal from the structure, 
on one hand, or induce further shell fragmentation, on the other 
hand. Many studies on this subject are present in the literature (see 
the corresponding discussion in ICiardi & Ferrar J20051 and refer- 
ences therein), but up to now comprehensive conclusions have not 
been reached: in fact the variety of approaches, approximations and 
parametrizations adopted by different authors leads to results which 
are often discordant. 



4.2 Clumping factor and recombination time 

Following Hernq uist & Springe! l2003h . we adopt a simple statis- 
tical description of the reionisation process. Noting that each ion- 
ising photon emitted is either absorbed by a newly ionised or by 
a recombining hydrogen atom, the filling factor Q(t) of regions 
of ionised hydrogen at any given time t follows from subtracting 
the total number of recombinations per atom from the total num- 
ber of ionising photons per hydrogen atom emitted before t. In this 
way, it is possible to derive a simple differential equation describing 



the transition from the neutral to the completely reionised universe 
ftVTadau et alJl999h . 



(6) 



c!Q _ n ioll _ Q_ 

dt UH ircc 

where uh represents the mean comoving density of hydrogen 
atoms and hion is the comoving emission rate of photons capable 
of ionising hydrogen. The recombination time of hydrogen, t rec , is 

1 



^rcc — 



(7) 



a B C n H {I + z)'i (I + 2 X ) ' 

where olb the recombination coefficient and \ is the relative abun- 
dance of helium with respect to hydrogen. The dimensionless 
clumping factor C of hydrogen is given by 



C 



{n 2 H+) 



(8) 



It is evident from Eq. that the recombination time is a function 
of redshift. Assuming a gas temperature of 10 4 K above which line 
cooling by atomic hydrogen becomes efficient, the previous rela- 
tion can be written as 



^rec (2) 



x 10 J 



1 



Myr 



(9) 



(i + zr 

(Mad auetalll999h . 

The left panel in Fig. [6] shows the redshift evolution of the 
clumping factor C, which has been estimated in our SPH simula- 
tions to be 



C 



Yji '»•/'. : Y. , "<;I>J 



(10) 



ISpringel & Hernquist! Eool . To avoid numerical artifacts pro- 
duced by gas particles belonging to formed objects which would 
artificially increase the value of C due to their high density, the 
sum extends over all gas particles whose density pi is smaller than 
a given threshold. In order to bracket realistic values for the thresh- 
old, we show the results (by shaded regions) assuming a maximum 
overdensity in the range between 100 and 500. We tested that the 
following analysis depends only slightly on the exact threshold. 

In general the growth of C with time is quite regular for all 
models. For this reason we decided to fit its behaviour assuming a 
parabolic relation, C = an + aiZ+ (I2Z 2 . The resulting best-fitting 
parameters 04, obtained using for C the central values of the shaded 
regions with the corresponding dispersion, are reported in TableQ 
Notice that the relation for the ACDM-RSI model is valid only for 
2 < 22, because of the strong flattening at higher redshifts. In line 
with our previous results, we find that gas clumping decreases from 
SUGRA through ACDM, SUGRA-RSI and ACDM-RSI. The dom- 
inant effect is caused by the shape of primordial power spectrum. 
In fact, the differences between the two RSI models are very small 
at all redshifts, where C is always between unity (corresponding 
to an almost perfectly smooth gas distribution), and three. For the 
ACDM model instead, the clumping factor starts to significantly 
deviate from unity at 2 ~ 25, reaching four at z ~ 18. As men- 
tioned before, the gas distribution appears much more clumpy in 
SUGRA, where C starts to grow already from z ~ 30, reaching 
C ~ 6 at 2 ~ 19. Notice that our values for C in ACDM are larger 
by 50 per cent compared to those obtained for the same model by 
Springel & Hernauist 1 2003 ) (see the open squares and the dashed 
line in the plot). This discrepancy is caused both by the higher res- 
olution of our simulation and by the different treatment of the gas 
component in the numerical code. Even if star formation is not 
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Figure 6. Left panel: The redshift evolution of the clumping factor C is shown for the different models (SUGRA, ACDM, SUGRA-RSI and ACDM-RSI 
from top to bottom). The shaded regions represent the uncertainty due to different choices for the maximum overdensity (here assumed between 100, dotted 
lines, and 500, solid lines; see the text for more detail). For reference, open squares and the dashed line show the results obtained for the Z4 model in 
ISpringel & H ernauist 1 2003). The right panel shows the evolution of the recombination time t ICC in Myr for the same models. 



Table 1. The fitting parameters for the redshift evolution of the clumping 
factor C computed in the different simulations (listed in Column 1). We 
adopt a parabolic relation, given by C = ao + a\z + a,2Z 2 . 



Model 


ao 


ai 


a 2 


SUGRA 


30.379 


-1.858 


0.030 


ACDM 


17.295 


-1.093 


0.019 


SUGRA-RSI 


9.592 


-0.585 


0.010 


ACDM-RSI 


20.086 


-1.655 


0.037 



included like in SDrinael & Hernauist 120031) . the fluid quantities 
(density, temperature and chemical abundances) are accurately and 
self-consistently evolved in our simulations, in particular also the 
cooling by molecular hydrogen. Thus, they are expected to produce 
a more realistic gas distribution in the diffuse intergalactic medium. 

Knowing the dumpiness factor C, it is easy to estimate the 
recombination time t Tec from Eq. |9j. The corresponding results 
are plotted in the right panel of Fig. |6| Of course, the maximum 
is obtained for the model with the lowest C, i.e. ACDM-RSI, in 
which the recombination process requires the long time of almost 
40 Myr at z ~ 18. The values for t ICC slightly decrease to about 
30 Myr for SUGRA-RSI and to 20 Myr for ACDM at the same 
redshift, and to less than 10 Myr at z ~ 19 for SUGRA model. The 
differences in both C and t xec are again spread over almost an order 
of magnitude because they are integrated observables like n c iouds- 



5 IMPLICATIONS FOR REIONISATION 

It is interesting to draw consequences from our results regarding 
the phenomenology of the overall reionisation process, checking in 
particular how these scenarios might be constrained by data on the 
reionised optical depth. Even if more accurate and quantitative esti- 
mates would require detailed simulations including radiative trans- 
fer, we can reliably discuss the problem of reionisation based on 
preceding work. 

We first need a relation between the reionisation epoch and 
the total Thomson optical depth r, which can be derived from the 
CMB data. To compare the WMAP results to the predictions for 
the cosmological model here considered, we compute the redshift 
evolution of r, adopting a simple model which assumes that com- 
plete reionisation occurs instantaneously at a some redshift z. The 
differential Thomson optical depth for complete ionisation is given 
by 

n(z) ot c 



dr{z) 
where 

ot 



(l + z)H(z) 



dz , 



Stt 
3 



6.65 x 10" 25 cm 2 



(ID 



(12) 



is the Thomson cross section, c is the speed of light and H(z) 
is Hubble parameter at z. The mean electron number density is 
given in terms of the ionisation fraction x D (z) by n(z) — n(z = 
0)(1 + z) A Xe{z). Notice that in the previous equations the only de- 
pendence on the cosmological model appears in H(z) through the 
Hubble constant and the expansion rates. For this reason the optical 
depth is completely independent on the spectral parameters (such 
as the power-spectrum index ri), and the RSI predictions agree ex- 
actly with the corresponding n = 1 model. 

The integrated optical depth, computed assuming x e = 1 
always, is shown in Fig. [7] for both ACDM and SUGRA with 
wo = —0.85. The horizontal dotted lines represent the observa- 
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Figure 7. The integrated Thompson optical depth r for ACDM (solid line) 
and SUGRA model with wq = —0.85 (long-dashed line) are compared to 
the WMAP constraints: first-year and three-years results are shown by the 
regions within the dotted lines and by the shaded region, respectively. 

tional estimat e (with its l-a error bars) obtained from the first-year 
WMAP data iKogut et alJl2003f) . while the shaded re gion shows 
the m ore recent values from the three-year WMAP data I Pa ge et all 
2006). 2 When compared to standard ACDM, SUGRA has a lower 
optical depth, with a difference of about At = 0.02 at z ~ 17 
and At = 0.01 at z ~ 10. As expected, this implies an ear- 
lier reionisation epoch, causing tension between the present model 
and the measured reionisation optical depth <Spergeletail E006). 
in addition to that regarding other cosmological probes (see, e.g., 
ISpergel et all2003tlKnop et all2003t|Riess et all2004h . Although 
a more accurate estimate would require a realistic description of 
the reionisation history, which in turn strongly depends on the free 
parameters of the adopted modelling, we can then assume that the 
reionisation in both models (ACDM and SUGRA) must occur and 
conclude at early epoch, ranging between z ~ 15 and z ~ 10, to 
agree with the measured optical depth within 1-cr. In general, these 
results demonstrate that the records of the high-redshift expansion 
rate in the physics of the primordial gas cloud s may be used to con- 
strain the dark energy from reionisation data. lMainini et all 120031) 
came to similar conclusions without considering modifications to 
the physics of gas clouds. 

We can now study in more detail how the reionisation pro- 
cess is going on in the redshift range covered by our simulations. 
First, we compute the total number of ionising photons nion. For 
this goal, we can follow the usual 'one star pe r halo' assumption 
for mini-halos (see, e.g.. lYoshida et al]l003albl) and safely use the 
cloud number as a good approximation for the number of very mas- 
sive stars: in fact strong radiative feedback disfavours the forma- 
tion of multiple stars inside each primordial gas cloud. Then we set 
the mass of a P op III star to be 500 Mq and we use the tables in 
ISchaered 12002) to estimate its lifetime (~ 1.9 Myr) and the mean 

2 We have to notice that our results for t must be considered as an extrapo- 
lation at redshifts higher than the reionisation epoch suggested by WMAP: 
since x e vanishes, the curves should indeed display a horizontal floor. 



ionising flux along the evolutionary track (~ 6.8 x 10 50 photons 
per second). We also assume an optimistic constant photon escape 
fraction of unity, which, however, can be considered realistic when 
the gas distributi on is reasonably smo oth lOh et all2001l) . Detailed 
calculations bv lKitavama et alJ l2004) indeed find that such a large 
escape fraction is plausible for small mass haloes. In order to com- 
plete the reionisation process, the resulting values for m on must be 
at least as large as the total number of hydrogen atoms in the simu- 
lation volume, which in our case corresponds to n to t ~ 4 x 10 66 . 
In Fig.[8|the dashed lines show the redshift evolution for rii on for 
the different models. We notice that only the SUGRA model is able 
to produce a sufficient number of ionising photons in the redshift 
range considered by our simulations, while for all remaining mod- 
els the ratio nj on /ntot is smaller than unity. In particular the two 
RSI models at the final simulation outputs reach only the 10-15 per 
cent of the required photons. 

Second, we have to consider that, as discussed before, recom- 
bination is counteracting the reionisation process. Using the values 
for the clumping factor and recombination time accurately derived 
from the simulations (and shown in Fig.|6j, we can easily compute 
the cumulative number of recombined atoms, numerically solving 
Eq. |6j. The results are also displayed in Fig. [3] where the shaded 
regions correspond to the uncertainties for C and t rec discussed 
in the previous section. The plot is suggesting that the recombina- 
tion process is rapid enough to significantly affect the fraction of 
ionised atoms. At the final simulation outputs, we find that the ra- 
tio between recombined and ionised atoms is about 75 per cent for 
SUGRA and ACDM, and about 25 for the two RSI models. 

Our analysis suggests that all models here considered, with the 
exception of SUGRA, are producing an insufficient number of col- 
lapsed objects to achieve complete reionisation in the redshift in- 
terval covered by the simulations. This is true not only for both RSI 
models, whose suppression of small scale power prevents the reion- 
isation, but also for the standard ACDM model. These conclusions 
ar e in qualitative agree ment with the previous analysis performed 
bv lYoshida et alj l003tl) . who, considering a simple model based 
on the formation of first stars in mini-haloes by molecular cooling, 
demonstrated that complete reionisation is reached when a density 
of 50 ~ 100 very massive stars per comoving Mpc 3 is turn ed on 
per ~ 1 recombination time (see also Yoshida et al. 003a cj). We 
notice that for these models the results are not conflicting with the 
three-year WMAP data which suggest a later reionisation epoch 
compared to the first-year analysis. Allowing the simulations to 
evolve up to z ~ 10 would certainly increase both the number of 
gas clouds and the dumpiness, completing the reionisation process. 

Finally, we notice that for the formation and evolution of cos- 
mic structures in SUGRA are so fast that the corresponding num- 
ber of gas clouds is so large and t ICC so low that a global reioni- 
sation is expected at redshifts even higher than z = 20. As men- 
tioned already, this anticipation of the reionisation epoch could be 
in confl ict with observations (see, e.g., the contour levels in Fig. 3 
of lSpergeletaiT2 006). indicating the potential constraining power 
of the dark energy records in the formation of primordial gas clouds 
which we pointed out in this work. 



6 CONCLUSIONS 

We presented hydrodynamic N-body simulations of the formation 
of primordial gas clouds on scales of tens of kpc in a variety of 
cosmological models, characterised by different dynamics in the 
dark-energy component. Our main results are that the records of 
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Figure 8. The redshift evolution of the number of ionizing photons n; on 
(dashed lines) and recombined atoms (shaded regions corresponding to the 
uncertainties for t rcc displayed in the right panel of Fig. |6| are shown for 
the different models (SUGRA, ACDM, SUGRA-RSI and ACDM-RSI from 
top to bottom). Both quantities are normalized to the total number of atoms 
ntot (see more detail in the text). 



the modified expansion rate are well evident in the population and 
the dumpiness of such clouds. Cosmological models with the same 
power-spectrum normalisation at present show earlier cloud forma- 
tion if the dynamics of the dark energy is enhanced, represented by 
an equation of state parameter w > — 1 as in quintessence models. 

Within dark energy models compatible with the present data 
on cosmic microwave background and large scale structure, the dif- 
ference in the integral population of clouds may vary by up to an 
order of magnitude, as a consequence of the different differential 
efficiency for structure formation. This is consistent with earlier re- 
sults indicating a h igher concentratio n in dark matter haloes under 
similar conditions jPolag et all2 004). 

Since abundance and dumpiness of structures are directly 
related to the amount of primordial power on the corresponding 
scales, we varied the shape of the primordial power spectrum by 
a running spectral index reducing power on small scales within 
the co nfidence level of the three-years WMAP data ISpergel et alJ 
2006). As expected, we find that the extra population and dumpi- 
ness of clouds produced by a higher dark energy abundance com- 
pared with its level today might be mitigated if the primordial spec- 
tral index is running, decreasing the power on small scales. 

Adopting a simple picture for the reionisation process, we de- 
rived consequences for the reionisation itself, leading to an earlier 
beginning of the reionisation process in models where cloud for- 
mation starts earlier. On the basis of these results, we are able to 
identify possible tension between the WMAP data on the reionisa- 
tion optical depth and cosmological models whose dark energy is 
as dynamical as in SUGRA quintessence models. 

Our results demonstrate that the effects on cosmological struc- 
ture formation from a modified expansion history through different 
dark energy models must be traced back to the formation of the first 
clouds. In turn, this means that constraints on the dark energy, and 
in particular its abundance at high redshifts, may be obtained by 



forthcoming experiments aiming at measuring the abundance and 
the dumpiness of primordial gas clouds. In particular, the Atacama 
Large Millimeter Array (ALMA 3 ) and the Mileura Widefield Ar- 
ray (MWA 4 ) will probe the early stages of structure formation, say 
between z = 6 and 10, where the reionization in progress should 
keep a record of the population of reionizing primordial gas clouds. 
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